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I. MODEL OF SUM OF LORENTZIANS 


We begin our analysis with fitting the optical data with a model dielectric function of a 
sum of Lorentzians. By htting experimental reflection and transmission spectra with this 
model, we can get separately temperature dependence of the spectral weight and scattering 
rate of the free carriers (Drude term), as well as resonance frequencies, spectral weight, 
and scattering rate for phonons and electronic transitions. For the present system, we 
use a Drude term, seven expected phonons, and hve electronic transitions. Five Lorentzian 
oscillators were found sufficient to give a very good representation of the electronic response. 
Figure Sl(a) shows the measured reflectivity and transmission spectra of Eu 2 lr 207 in the 
frequency range of phonons together with the model curves. These are the same spectra as 
in hgure 1 of the paper but with linear frequency scale. 


Figure Sl(b) shows each Lorentzian term separately except for the phonons which are 
shown as one spectrum. The sum of these terms gives the net optical conductivity. It is 
clear that these components of the optical conductivity have different resonance frequencies 
and,especially, widths and that they do not interfere with each other in the htting procedure. 
We note that the Drude term at 7 K (blue curve) is broad and it gives a small but hnite 
dc conductivity. The fact that fact that Drude term (red curve) and the lowest electronic 
peak (magenta curve) have opposite slope allows the clean separation of their contributions 
by htting the transmission spectra. We hrst ht simultaneously the 7 K transmission and 
rehection spectra. At this temperature, the transmission is maximum at low frequencies 
and it is a clear signature of semimetallic optical conductivity. We found that allowing for a 
temperature dependence of the Drude and phonon parameters is sufficient to account for the 
observed temperature dependence. These ht results are shown in Fig. 3 of the paper as open 
circles. We then ht the rehectivity spectra only at temperatures at 70 K and higher where 
there is no transmission for the thinned crstal at these temperatures. Again, a variation 
of only the Drude and phonon parameters keeping all other parameters hxed is sufficient. 
Results of this ht are shown in Fig. 3 as solid squares. 
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FIG. 1. (color online), (a) Reflectivity (R) and transmission (Tra) of a Eu 2 lr 207 single crystal. 
Dashed curves are fits using the Lorentzian model (Eq. (1) of the paper), (b) Components of the 
total optical conductivity: Drude, phonons, and three lowest in frequency electronic transitions. 


II. MODEL OF THE VARIATIONAL DIELECTRIC FUNCTION (VDF) 


In the second method, we use the results of the hrst method as a starting spectrum and 
calculate the VDF ^(ci;) [27]. This accounts for small deviations between Loretzian £t and 
experimental data. Figure S2 shows an example of the two types of £t. It is apparent 
from the hgure that the model curve (blue) for VDF £t is practically indistinguishable from 
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FIG. 2. (color online). Experimental and model reflectivity of a Eu 2 lr 207 single crystal calculated 
from the full variational dielectric function. 


the experimental curve (red). The VDF ht is equivalent to the standard Kramers-Kronig 
transformation technique [27] and both rely on assumed extrapolation beyond the measured 
interval. The value of the VDF £t is that it reproduces all small features of the spectrum 
not included in the sum of Lorentzians. 


III. SPECTRAL WEIGHT 

To check our conclusion that all temperature changes come only from the Drude and 
phonon terms, we use another standard spectroscopic tool — the analysis of the total spectral 
weight SW{uj)\ 

SW{uc) = ai{uj)du = (1) 

where ai is the optical conductivity in units of (Ohm-cm)“^ and the effective (or total) 
plasma frequency ojpt includes all possible excitations. In Figure S3, we compare obtained 
by integration of the VDF optical conductivity up to the cutoff frequency a;c=530 cm”^ 
with the Drude weight, where the oJsd data are the same as in Fig. 3(a) of the paper 
obtained from the sum of Lorentzians £t. Here, we have shifted the curve (blue stars) 
by an arbitrary constant to exclude all temperature independent excitations. We hnd that 
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FIG. 3. (color online). Comparison of the temperature dependences of the total (stars) and Drude 
(open circles and solid squares) spectral weights. The total spectral weight was shifted down to 
exclude the temperature independent parts. 

the agreement of all data is quite impressive considering that not only are they calculated 
differently but also they are based on different measured spectra: blue stars were obtained 
from the reflectivity measured with the 4 K bolometer in the frequency range 90-600 cm“^, 
open circles and solid squares come from the transmission and reflection measured with the 
1.6 K bolometer in the frequency range 10-200 cm“^. From a comparison of the u}p^{T) (total 
spectral weight) and the a;^^(T) (Drude term), we conclude that only the Drude spectral 
weight depends on temperature. While the phonons do have temperature dependent widths 
their total spectral weight is practically independent of temperature compared to the Drude 
weight. 
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A Weyl semimetallic state with pairs of nondegenerate Dirac cones in three dimensions was 
recently predicted to occur in the antiferromagnetic state of the pyrochlore iridates. Here, we 
show that the THz optical conductivity and temperature dependence of the free carrier response 
in pyrochlore Eu 2 lr 207 match the predictions for a Weyl semimetal and suggest novel Dirac liquid 
behavior. The interband optical conductivity vanishes continuously at low frequencies signifying 
a semimetal. The metal-semimetal transition at Tn = 110 K is manifested in the Drude spectral 
weight, which is independent of temperature in the metallic phase, and which decreases smoothly 
in the ordered phase. The temperature dependence of the free carrier weight below Tat is in good 
agreement with theoretical predictions for a Dirac material. The data yield a Fermi velocity vp ks 
4 ■ 10^ cm/s, a logarithmic renormalization scale An ~ 600 K, and require a Fermi temperature 
of Tp « 100 K associated with residual unintentional doping to account for the low temperature 
optical response and dc resistivity. 

PACS numbers: 75.47.Lx, 76.30.Kg, 78.30.-j, 71.30.-|-h 


Condensed matter physics is in an exciting new era 
dominated by Dirac materials and topological effects. 
The introduction of topology, particularly in association 
with spin-orbit coupling, has led to predictions of novel 
electronic and optical phenomenai^^— . The most recent 
development is the discovery of three-dimensional ( 3 D) 
Dirao^^— and Weyliip— materials, where the Dirac cones 
come in pairs and have non-Kramers degenerate chiral 
bands protected by topology. A 3 D Weyl state has been 
predicted^^ for the 5 d transition metal oxide family of 
pyrochlore iridatesi^, which are strongly interacting ma¬ 
terials with a strong spin-orbit interaction. However, in 
the absence of direct experimental evidence their ground 
state is still under intense discussioiJ^^— . 

The pyrochlore iridates, i?2li'207, are strongly in¬ 
teracting materials that exhibit frustrated magnetism 
and an associated metal-insulator transition as R varies 
across the rare earth series^S. Due to a strong spin- 
orbit interaction these materials have been recognized 
as having potentially exotic ground states including the 
Weyl semimetallic state in the low temperature magnetic 
phase. Also, axion insulato r^"^d^ , topological band in¬ 
sulator and Mott insulator—, and spin liquid states^ii^ 
have been proposed. X-ray diffraction experiments show 
that there is no structural transition breaking the cubic 
symmetry in Nd, Eu, and Pr pyrochlores, and that Eu 
compound displays a smooth thermal contraction of the 
lattice parameters through the ” metal-insulator” transi- 
tion^^. Ueda et ah— reported a 50 meV optical gap in 
the ground state of polycrystalline Nd2lr207 which was 
tuned to zero by the partial substitution of Ir by Rh and 
that gapless state was proposed to be a Weyl semimetal. 
By contrast, infrared studies of Bi2lr207 found a strongly 


metallic ground state^^. Therefore, at present there is no 
clear evidence of the predicted Weyl semimetal state in 
the rare earth pyrochlores. 

Here, we present evidence of a Weyl semimetal on the 
basis of experimental and theoretical studies of the opti¬ 
cal response, and we identify the electronic phase transi¬ 
tion in Eu2lr207 at T/v as a metal-semimetal transition in 
which the conductance is controled by the thermal pop¬ 
ulation of the Weyl cones. Single crystals of Eu2lr207 
were grown as described elsewhere^^. Our crystal has 
the same dc conductivity as the sample #2 of Refj^. 
Fourier transform infrared reflectivity and transmission 
measurements were performed on the 1.8 mm in diameter 
( 111 ) face of a single crystal. For transmission measure¬ 
ments at frequencies below 100 cm“^, the crystal was 
glued down to an intrinsic Si substrate using transparent 
Stycast epoxy and was polished down to ~ 10 micron. 

Figure 1 shows the measured reflectivity and transmis¬ 
sion spectra of Eu2lr207 over a broad frequency range. 
The transmission measured at 7 K immediately suggests 
a semimetalic state. The negative slope implies an op¬ 
tical conductivity growing with frequency. An insulator 
(or gapped state) would give a flat transmission, and a 
free carrier Drude term would give a slope of opposite 
sign. 

We have analyzed the reflectivity and transmission 
spectra in two standard ways^l: first, using a model 
of a sum of Lorentzian oscillators and second, using a 
variational dielectric function (equivalent to a Kramers- 
Kronig transformation)^. In the first method, we fit the 
reflectivity and the transmission spectra of Eu2lr207 us¬ 
ing a Lorentzian model, in which the complex dielectric 
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FIG. 1. (color online). Reflectivity (R) and transmission 
(Tra) of a Eu 2 lr 207 single crystal. Dashed curves are fits 
using the Lorentzian model (Eq. ([T])). Grey lines are spectra 
at intermediate temperatures. The temperature dependence 
was measured over a 20-6,000 cm“^ frequency range. 


function e = ei + ie 2 takes the form: 
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e{uj) = 




• lUJ'Jj 


( 1 ) 


Here woj, 7j, and ujgj are the resonance, the damping, 
and the spectral weight frequencies of the jth electric 
dipole active mode, respectively. For the present system, 
we use a Drude term, seven symmetry allowed phonons^S, 
and we find that the electronic response can be well repre¬ 
sented using five electronic transitions, the lowest one at 
a;o=306 cm~^ with 7=1143 cm“^. In the second method, 
we use the results of the first method as a starting spec¬ 
trum for the variational e(w) to account for small devi- 
ations^i. The complex optical conductivity is related to 
the dielectric function by cr = cti -I- ia 2 = uje/Airi. 

Figure 2 presents the frequency-dependent optical 
conductivity ui and the dielectric constant Si as ob¬ 
tained from the analysis of the experimental spectra. 
Figure 2(a) shows the optical conductivity in a broad fre¬ 
quency range as obtained from the variational e analysis. 
The two-peak conductivity spectra are consistent with 
resonant inelastic X-ray scattering spectra^ and band 
structure calculationsi^. The vanishing of cti at low fre¬ 
quencies is a signature of a Dirac semimetal. The 1 eV 
band consists of transitions within the Jeff=1/2 bands^^. 
Figures 2(b) and 2(c) show the optical conductivity and 
the dielectric constant at low frequencies. The increasing 
£i at low frequencies (Fig. 2(c)) is another characteristic 
of a semimetal and is expected to diverge logarithmically 
for an intrinsic Weyl semimetal in the absence of dis¬ 
order—. The optical conductivity of an intrinsic Weyl 
semimetal can be written as 


A . ^ 

12 h vf 


( 2 ) 


FIG. 2. (color online). Optical conductivity and dielec¬ 
tric constant of a Eu 2 lr 207 single crystal, (a) Temperature- 
independent broad band electronic optical conductivity, (b) 
Temperature-dependent low-frequency part of the optical 
conductivity. Our zero frequency values are close to the dc 
measurements^^. The magenta line is a linear fit of equa¬ 
tion m to the conductivity, (c) The dielectric constant. The 
continuous green line is the blue 7 K line with the phonon 
contributions removed. The dashed black line denotes the 
high-frequency dielectric constant k obtained from the plas- 
mon fit. 


arising from interband excitations within the Dirac 
conei^. Here, we denote by g the degeneracy of Dirac 
nodes and by vp = {vpVpVpy^^ the geometric mean of 
the Fermi velocities. By fitting this expression to the 
optical conductivity (magenta line in Fig. 2(b)) using 
g=24ii, we obtain vp = 3.4 • 10^ cm/s. 

Based on this fit to the optical conductivity, we can ex¬ 
tract the Coulomb interaction strength in the semimetal 
as given by the effective Dirac fine structure constant 
a=e^/hvpK. Here, k which is of the order of 10 is 
the high-frequency dielectric constant which shows no 
temperature-dependence as is evident in Fig. 2(c). (We 
use e and k to designate experimental and theoretical di¬ 
electric constant, respectively.) This gives a bare value of 
0=0.7 for the interaction strength, indicating that elec¬ 
tron interaction could have considerable impact on the 
low-energy properties of this Dirac material. Since the 
corresponding fine structure constant for quantum elec¬ 
trodynamics is 1/137, the current problem is a strong¬ 
coupling (and non-relativistic) version of quantum elec¬ 
trodynamics. 

The conductivity exhibits a temperature dependence 
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only at low frequencies ti;<500 cm Fits of the Eq. (HD 
to the experimental spectra reveal that only the Drude 
and phonon terms are temperature dependent while the 
interband conductivity can be kept constant at all tem¬ 
peratures. The validity of this conclusion is verified by 
a spectral weight analysis presented in the Supplemen¬ 
tal Material^l. We show the result for the experimental 
Drude spectral weight and scattering 7 as a function of 
temperature in Fig. 3(a,b). The metal-semimetal tran¬ 
sition at T/v is manifested as a smooth decrease of the 
Drude spectral weight in the Weyl state which can exist 
only below where time- reversal symmetry is bro- 
keni^. 

In Fig. 3(a), we compare the experimental Drude 
weight with the theoretical prediction for a Dirac liq¬ 
uid. The low-temperature saturation indicates a finite 
Fermi energy Ep (associated with unintentional resid¬ 
ual doping), leads to a low temperature dc conductiv¬ 
ity and marks a transition between an extrinsic metallic 
low-temperature regime {ksT < Ep) and an intrinsic 
semimetallic high-temperature regime {ksT > Ep). Our 
observations are consistent with the Drude weight as pre¬ 
dicted by kinetic theory 


2,-i2 


e V 


^sD — 


dED{E) 


df{E) 
dE '■ 


(3) 


where D{E) = gE'^/[2Tr‘^{hvp)^] is the density of states 
and f{e) the Fermi-Dirac distribution. A fit of the 
Drude weight (solid black line in Fig. 3(a)) gives a fi¬ 
nite Fermi energy of Ep/kp^lO K and a Fermi velocity 
of vp=4:-10'^ cm/s, in agreement with the direct fit to 
the frequency-dependent conductivity. For comparison, 
we include a fit to the high-temperature Drude weight in 
the strictly intrinsic limit (dashed blue line in Fig. 3(a)) 
which shows that finite doping does not affect the Drude 
weight at high temperature {kpT > Ep). Note that for 
a general semimetallic dispersion E ~ \k\^, the Drude 
weight is expected to scale as uj^p, ^ j’H-i/z tem- 
perature^^. The superlinear temperature dependence im¬ 
plies z < 1 in this effective model, compellingly exclud¬ 
ing, for example, a parabolic semimetallic ground state. 

Our experimental Drude weight is seen to deviate from 
the strict linear temperature scaling of a noninteract¬ 
ing semimetal, which can be attributed to electron in¬ 
teractions: the interaction strength of an intrinsic Dirac 
semimetals acquires a scale-dependence due to ultravi¬ 
olet renormalization, giving rise to a superlinear tem¬ 
perature dependence of the Drude weight and provid¬ 
ing a direct signature of electron interaction effect a^^'^^ . 
Within RPA, the temperature dependence of the inter¬ 
action strength is given by^i^i 

a(T) = -ln-i^. (4) 

9 T 

Instead of depending on e^, up-, and k separately, this ex¬ 
pression contains a single renormalized energy scale Ap, 
the Landau pole, at which the coupling a diverges. This 



FIG. 3. (color online). Temperature dependence of the spec¬ 
tral Drude weight and scattering rate 7 . (a) Open circles 

— results of simultaneous fit of reflectivity and transmission, 
solid squares — results of reflectivity fits, both using Eq. 0 , 
lines — fits of kinetic theory (Eq. ([3))) and plasmon to these 
experimental points below Tjv. (b) Experimental scattering 
rate of free carriers. 


Landau scale is an effective parameter of the system that 
can be extracted from measurements of the Drude weight 
and corresponds to a cutoff scale beyond which the Dirac 
dispersion is no longer linear. In the present system, it 
can be associated with the Lifshitz saddle point between 
each pair of Weyl points^. 

In the absence of disorder and phonon contributions, 
the Drude weight corresponds to the plasmon frequency 
as defined by the zero of the dielectric function e{uj) = 
Ko — . We assume kq independent of T and uj 

(Fig. 2(c)) and it is a good assumption in the frequency 
range 300-800 cm“^ where solutions for lOsd are found. 
We perform a fit to the experimental uJsd using the plas¬ 
mon dispersion of an extrinsic Dirac semimetal^, which 
fully takes into account the competing effects of finite 
doping and renormalization of the interaction. Since the¬ 
oretical calculations only predict the renormalization of 
the effective fine structure constant a, there is a question 
whether renormalization effects should be manifested in 
the Drude weight, the dielectric constant, or both. The 
result of the fit is shown in Fig. 3(a) as a red solid line. 
We obtain good agreement with the experiment with an 
effective static dielectric constant of ko« 13 and a Fermi 
energy of Ep/kp^WQ K, which is consistent with the 
experimental dielectric function ei (Fig. 2), and the fit 
to the Drude-Boltzmann weight (Eq. [3]), respectively. 
The Landau pole takes an anomalously small value of 
AiSibOO K, giving rise to strong superlinear scaling. This 
could be taken as an indirect sign that the material hosts 
pairs of Weyl cones in close proximity, consistent with 
the picture of a Weyl semimetal state with broken time- 
reversal symmetry. For comparison, we also include a fit 
of the intrinsic finite-temperature plasmon in Fig. 3(a) 
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(dash-dotted magenta line). As before, the intrinsic re¬ 
sult captures the high-temperature behavior where the 
Weyl state is most expected and the saturation at low 
temperature can be attributed to finite doping. Intu¬ 
itively, one could expect a characteristic step in ai{ui) at 
2 A_fss 140 cm“^ but it should be smeared out by large 
scattering rate and cannot be seen clearly in Fig. 2(b). 

Linear (or quasi-linear) dependence of the low fre¬ 
quency optical conductivity — a signature of a Dirac 
semimetallic state — has also been recently reported in 
several other 3D materials: HgCdTe^, quasicrystals^, 
and ZrTeg^. A notable result of our work compared 
to previous studies is the excellent agreement of the ex¬ 
perimental and theoretical temperature dependent Drude 
spectral weight providing strong evidence of a Weyl state 
with finite doping. However, the independence of the 
Drude parameters on temperature in the paramagnetic 
state above T^r is surprising and merits further study. 

In summary, the optical response of a Eu 2 lr 207 single 
crystal reveals a semimetallic electronic structure with 
approximately linear frequency dependence of the op¬ 
tical conductivity down to 3 meV at low temperature. 
Below T/v, the Drude spectral weight diminishes con¬ 
sistent with the reduced thermal excitations of a Weyl 
semimetal. This means that the long-thought “metal- 
insulator” transition in pyrochlore iridates may actually 


be a metal-semimetal transition, at least, in the case of 
Eu 2 lr 207 . These two data sets can be modeled assum¬ 
ing a Weyl state, as 24 Weyl cones with an average Fermi 
velocity vf=4:-10'^ cm/s. The theoretical analyses of our 
optical data point toward signatures of the ultraviolet 
renormalization expected for an interacting Dirac liquid 
manifesting in the super-linear temperature dependence 
of the Drude weight. 
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